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Abstract The complete mitochondrial NADH dehydrogenase subunit 2 (ND2) gene sequences of two species of 
Thermophis, T. baileyi and T. zhaoermii, were sequenced and compared to those of 86 sequences from other snakes 
(74 from Caenophidia and 12 from Henophidia). By using Bayesian inference (BI) and maximum likehood (ML) 
approaches, Thermophis was demonstrated as the sister group to the North American relicts of Dipsadinae, and rooted in 
Central and South American members of this subfamily. The results suggest that the closest relatives of Thermophis are 
the North American relicts, and thus support the hypothesis for an Asian-North American origin of xenodontine snakes 
as suggested by Vidal et al. (2000). Extensive sampling of Asian snakes and American dipsadines is needed to further 


test this hypothesis in the future. 
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1. Introduction 


Hot-spring snakes are endemic to the Tibetan Plateau, 
known as “the roof of the world” and “the world's third 
pole”. Their precise systematic position is still a question. 

Initially, the hot-spring snakes such as Tropidonofus 
baileyi (= Natrix baileyi) from the Tibetan Plateau 
were placed within Natricinae (Wall, 1907). However, 
Malnate (1953) designated a new (monotypic) genus 
Thermophis within Xenodontinae for these snakes 
based on the characters from the hypapophysis and the 
hemipenis. Additionally, he considered Thermophis to 
closely resemble the ancestor of both the Natricinae 
and Colubrinae. Thermophis was placed in Colubridae 
incertae sedis on account of the hemipenial character of 
a divided sulcus in a single calyculate lobe (described by 
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Malnate, 1953) (Zaher, 1999). The molecular study of 
Lawson et al. (2005) retained this classification based on 
the previous studies. 

Interspecific taxonomy within this genus and estimates 
of range did not change for nearly 100 years. Thermophis 
baileyi was thought to occur only in Lhasa and the 
Gyantze region (Malnate, 1953). However, T. baileyi also 
occurred in the middle reaches of the Yarlung Zangbo 
River such as in Lhatse and Chushur (Rao, 2000). 
Furthermore, 13 locations were recorded for discovering 
T. baileyi in the central part of the Tibet Plateau (Dorge et 
al., 2007). Thermophis was also found in Litang County, 
Sichuan Province (Liu and Zhao, 2004). Most recently, 
after detailed morphological studies, the suggested 
populations from Litang were considered a new species 
and named as T. zhaoermii (Guo et al., 2008). Huang et 
al. (2009) designated the hot-spring snakes (endemic to 
the Tibetan Plateau) as a "relict" species, according to 
Pinou et al. (2004). This designation was due in part to 
their unique ecological adaptations to the environment 
and climate of the Tibetan Plateau, but also to the lack of 
closely related taxa in Asia. 

Using unique characters (plesiomorphies and/or 
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apomorphies) to characterize and group some snake taxa 
is difficult due to the limited morphological information 
that can be obtained due to their specialized and 
therefore simple body structure (Huang et al., 2009). 
Thus, molecular data appear to play a more important 
role in snake systematics. Over the last two decades, 
fragment(s) of mitochondrial DNA (mtDNA) and 
single copy of nuclear DNA have frequently been used 
as a versatile and effective marker to investigate the 
phylogenetic relationships among snakes (e. g., Kelly et 
al., 2003; Lawson et al., 2005; Burbrink and Lawson, 
2007; Huang et al., 2009; Vidal et al., 2010; Pyron et 
al., 2011). For example, He et al. (2009) used cyt b and 
c-mos, and Huang et al. (2009) selected cyt b and ND4 
to examine the systematic position of Thermophis within 
Caenophidia, and found that Thermophis appeared to be a 
member of the subfamily Dispadinae (then Xenodontinae; 
see Pyron et al., 2011). But their data failed to clarify 
the relationship between the New World xenodontines 
and Thermophis. Therefore, our goal of this paper is to 
provide sufficient data for a better understanding of the 
phylogenetic position of Thermophis among Colubroidea. 

In this study, we sequenced 3 Thermophis and used 
73 Caenophidian ND2 sequences from GenBank, as well 
as 12 henophiids as outgroups in a phylogenetic analysis 
with an attempt to clarify the phylogenetic relationship 
between Dipsadinae and Thermophis. The systematic 
nomenclature within Colubridae used in this paper was 
primarily based on the classification recommended 
by Lawson et al. (2005), which was derived from 
their phylogenetic analyses of cyt b and c-mos genes, 
as well as on the evidences obtained from previous 
work using morphological, immunological, and DNA 
sequence data. However, we simultaneously followed the 
recommendation by Pyron et al. (2011) in using the name 
Dipsadinae to refer to Xenodontinae (South American 
xenodontines + Central America dipsadines + North 
American relict snakes) due to their controversial status at 
the subfamilial level (e. g., Vidal et al., 2000; Lawson et 
al., 2005; Zaher et al., 2009). Much of the classification 
proposed by Lawson et al. (2005) has been adopted and 
corroborated by many contemporary workers in current 
literature (e. g., He et al., 2009; Huang et al., 2009; Ling 
et al., 2010). 


2. Materials and Methods 


2.1 Downloaded sequences and sampling information 


We sampled 87 species (88 sequences) representing mostly 
Alethinophidia (primitive + advanced snakes; Table 1). Of 
these, 86 sequences were retrieved from GenBank, and 2 
were unique sequences from 8 individuals of Thermophis: 
6 T. baileyi from Yangbajain (HS2008 18-2; N: 30.0638; 
E: 90.4876; 4100 m.), Xizang (Tibet), China and 2 T. 
zhaoermii form Litang (HS200818-1; N: 29.6921, E: 
100.3842; 3700m), Sichuan, China. The snakes collected 
for this study were humanely euthanized using lethal 
injection. The fresh liver tissues or tail tip biopsies were 
removed and immediately preserved in 95% ethanol for 
sequencing their mitochondrial ND2 gene. Specimens 
are deposited in the museums of Sichuan University and 
Sichuan Forestry Academy. 


2.2 DNA extraction, amplification and sequencing 
Eight individuals of Thermophis were collected during 
the summer of 2009. Three distal ends of their ventral 
scales were removed and immediately preserved in 
95% ethanol for sequencing their mitochondrial ND2. 
We extracted total genomic DNA using the phenol/ 
chloroform extraction procedure in Sambrook et al. 
(1989). Polymerase chain reaction (PCR) was used for 
amplification of mitochondrial ND2 with primer pairs 
L4437c, tRNA-trpR’ (Huang et al., 2007). PCR was 
carried out under the following conditions: an initial cycle 
at 94° for 7 min followed by 40 cycles of 94° denaturation 
(40 s each), 56° annealing (30 s) and 72° extension (1 
min) with a final extension at 72° for 7 min. 


2.3 Phylogenetic inference Newly obtained sequences 
were aligned with other retrieved sequences using 
Clustal X (Thomopson et al., 1997) and were translated 
into amino acids to check for the presence of stop 
codons within the sequences which might indicate the 
amplification of pseudogene(s). Nucleotide composition 
and the uncorrected P-distances were calculated in Mega 
4.0 (Tamura et al., 2007). 

Bayesian inference (BI) approaches were used to 
reconstruct phylogenetic bifurcating trees in the present 
study, using MrBayes V3.0 (Huelsenbeck and Ronquist, 
2001). Modeltest 3.7 (Posada and Crandall, 1998) was 
used to find the model of nucleotide evolution using the 
Bayesian Information Criterion (BIC) for BI, with the 
scores formed by PAUP* V4b10 (Swofford, 2002). In the 
Bayesian analysis, a starting tree was obtained using the 
neighbor-joining algorithm. Posterior probabilities (PP) 
were obtained by Markov Chain Monte Carlo (MCMC) 
analysis with one cold chain and three heated chains. 
Bayesian tree estimation and support for the nodes was 
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Table 1 Information of the sequences retrieved from GenBank and sequenced in this study. Much of the classification is proposed by Lawson 
et al. (2005) but simultaneously follows the recommendation by Pyron et al. (2011) in using the name Dipsadinae to refer to Xenodontinae 
(South American xenodontines + Central America dipsadines + North American relict snakes). 


Clade Genus and species Accession No. Source 
Colubridae 

Colubrinae Ptyas korros AY487023 Nagy et al. (2004) 
Bogertophis subocularis DQ902206 Burbrink and Lawson (Unpubl.) 
Pituophis deppei FJ627848 Pyron and Burbrink (2009) 
Cemophora coccinea DQ902249 Burbrink and Lawson (Unpubl.) 
Lampropeltis ruthveni FJ627847 Pyron and Burbrink (2009) 
Stilosoma extenuatum DQ902245 Burbrink and Lawson (Unpubl.) 
Pseudelaphe flavirufa DQ902216 Burbrink and Lawson (Unpubl.) 
Arizona elegans DQ902204 Burbrink and Lawson (Unpubl.) 
Pantherophis obsoletus FJ627844 Pyron and Burbrink (2009) 
Rhinocheilus lecontei FJ627838 Pyron and Burbrink (2009) 
Senticolis triaspis DQ902237 Burbrink and Lawson (Unpubl.) 
Coelognathus radiatus DQ902230 Burbrink and Lawson (Unpubl.) 
Oreophis porphyraceus DQ902226 Burbrink and Lawson (Unpubl.) 
Elaphe quatuorlineata AY487028 Nagy et al. (2004) 
Euprepiophis mandarinus DQ902222 Burbrink and Lawson (Unpubl.) 
Dinodon semicarinatus NC_001945 Kumazawa et al. (1998) 
Coronella girondica AY487027 Nagy et al. (2004) 
Gonyosoma oxycephalum DQ902241 Burbrink and Lawson (Unpubl.) 
Drymarchon corais DQ902207 Burbrink and Lawson (Unpubl.) 
Salvadora mexicana AY487036 Nagy et al. (2004) 
Masticophis flagellum piceus AY487022 Nagy et al. (2004) 
Coluber constrictor mormon AY487002 Nagy et al. (2004) 
Macroprotodon brevis DQ902270 Busack and Lawson (2008) 
Hemerophis socotrae AY487016 Nagy et al. (2004) 
Eirenis levantinus AY487032 Nagy et al. (2004) 
Hierophis viridiflavus AY487018 Nagy et al. (2004) 
Spalerosophis diadema AY487020 Nagy et al. (2004) 
Hemorrhois ravergieri AY487011 Nagy et al. (2004) 
Platyceps najadum AY487009 Nagy et al. (2004) 
Oocatochus rufodorsatus DQ902232 Burbrink and Lawson (Unpubl.) 
Orthriophis cantoris DQ902246 Burbrink and Lawson (Unpubl.) 

Xenodermatinae Achalinus meiguensis NC_011576 Wang et al. (Unpubl.) 

Homalopsinae Enhydris plumbea NC_010200 Yan et al. (2008) 

Natricinae Amphiesma sauteri AF384824 Alfaro and Arnold (2001) 
Natrix tessellata AY487776 Guicking et al. (2006) 
Virginia striatula AF384852 Alfaro and Arnold (2001) 
Tropidoclonion lineatum AF384850 Alfaro and Arnold (2001) 
Nerodia harteri AF384854 Alfaro and Arnold (2001) 
Xenochrophis flavipunctatus FJ416748 Hedges et al. (2009) 
Regina septemvittata AF384836 Alfaro and Arnold (2001) 
Clonophis kirtlandii AF384827 Alfaro and Arnold (2001) 
Thamnophis validus validus EF417447 de Queiroz and Lawson (2008) 
Seminatrix pygeae AF384839 Alfaro and Arnold (2001) 
Storeria occipitomaculata AF384840 Alfaro and Arnold (2001) 

Dipsadinae Hypsiglena EU728580 Mulcahy and Macey (2009) 
Pseudoleptodeira latifasciata EU728579 Mulcahy and Macey (2009) 
Imantodes cenchoa EU728586 Mulcahy and Macey (2009) 
Leptodeira septentrionalis EU728590 Mulcahy and Macey (2009) 
Sibon nebulatus EU728583 Mulcahy and Macey (2009) 
Helicops angulatus FJ416751 Hedges et al. (2009) 
Uromacer oxyrhynchus FJ416754 Hedges et al. (2009) 
Darlingtonia haetiana FJ416774 Hedges et al. (2009) 
Antillophis andreae FJ416781 Hedges et al. (2009) 
Alsophis cantherigerus ruttyi FJ416785 Hedges et al. (2009) 
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(Continued Table 1) 


Jaltris dorsalis FJ416773 Hedges et al. (2009) 
Arrhyton funereum FJ416777 Hedges et al. (2009) 
Hypsirhynchus ferox FJ416780 Hedges et al. (2009) 
Heterodon platirhinos FJ416750 Hedges et al. (2009) 
Farancia abacura DQ902239 Burbrink and Lawson (Unpubl.) 
Dipsadinae 
(incertae sedis) 
Thermophis baileyi JF411072 This study 
Thermophis zhaoermii JF411071 This study 
Thermophis zhaoermii NC_012816 He et al. (2009.) 
Acrochordidae Acrochordus granulatus NC_007400 Dong and Kumazawa (2005) 
Viperidae 
Viperinae Causus defilippi GU045452 Castoe et al. (2009a) 
Vipera aspis aspis AY581266 Garrigues et al. (2004; Unpubl.) 
Daboia russellii NC_ 011391 Chen and Fu (Unpubl.) 
Crotalinae Deinagkistrodon acutus DQ836185 Huang et al. (2007) 
Gloydius brevicaudus DQ839399 Huang et al. (2007) 
Crotalus viridis viridis AY704799 Wüster et al. (2005) 
Ovophis okinavensis AB175670 Dong and Kumazawa (2005) 
Agkistrodon piscivorus NC_009768 Jiang et al. (2007) 
Viridovipera stejnegeri stejnegeri FJ752492 Lin et al. (Unpubl.) 
Elapidae 
Bungarus multicinctus NC_011392 Chen and Fu (Unpubl.) 
Micrurus fulvius GU045453 Castoe et al. (2009b) 
Naja naja NC_010225 Yan et al. (2008) 
Ophiophagus hannah NC_011394 Chen and Fu (Unpubl.) 
Pythonidae 
Python regius NC_007399 Dong and Kumazawa (2005) 
Boidae 
Erycinae Charina trivirgata GQ200595 Douglas and Gower (2010) 
Boinae Eunectes notaeus AM236347 Douglas et al. (2006) 
Boa constrictor NC_007398 Dong and Kumazawa (2005) 
Xenopeltidae Xenopeltis unicolor NC_007402 Dong and Kumazawa (2005) 
Cylindrophiidae Cylindrophis ruffus AY662538 Townsendb et al. (2004) 
Uropeltidae Rhinophis philippious GQ200594 Douglas and Gower (2010) 
Tropidophiidae Tropidophis haetianus FJ755181 Castoe et al. (2009b) 
Aniliidae Anilius scytale GQ200593 Douglas and Gower (2010) 
Leptotyphlopidae Rena humilis AB079597 Kumazawa and Nishida (1995) 
Typhlopidae Typhlops reticulatus NC_010971 Castoe et al. (2008) 
Ramphotyphlops AM236346 Douglas et al. (2006) 
australis 


inferred using a Monte Carlo Markov chain model as 
implemented in MrBayes, with 10,000,000 generations 
after a "burn-in" of 3,000,000 generations. Sampling 
was performed every 1000 generations. We ran three 
additional analyses starting with random trees. Consensus 
of all the post-burn generations (3,600,000 generations 
resulting in 36,000 trees) was computed from all four 
runs. In addition, ML tree was derived using a RAxML 
BlackBox in the model of Maximum likelihood search 
(http:// phylobench.vital-it.ch/raxml-bb). 


3. Result 


3.1 Sequence characteristics The complete ND2 gene 


sequences (1032 bp, i. e., 344 codons) of T. baileyi and 
T. zhaoermii were obtained in this study, both including 
the stop codon TAG located at the 3 terminal end. An 
open reading frame was obtained for all ND2 sequences. 
The average nucleotide composition (%), numbers 
and percentages of variable (V) and PI sites for ND2 
gene fragments of Caenophidia, Dipsadinae (excluding 
Thermophis, clustered in this subfamily, see below), 
and Thermophis used in the present study were listed in 
Table 2. The percentages of the four nucleotides were not 
significantly different across taxa groups, and showed 
strong bias against guanidine, typical for mitochondrial 
DNA (e. g., Lawson et al., 2005; Huang et al., 2009). 

The uncorrected P-distance within Caenophidia, 
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Table 2 Average nucleotide composition (%), numbers and percentages (%) of variable (V) and parsimony informative (PI) sites for ND2 


gene, based on the data used in the present study. 


Taxa T(U) C A 

Caenophidia n=73 22.8 30.5 37.8 
Dipsadinae™ n=15 22.7 31.1 37.3 
Thermophis n=3 21.0 33.5 36.4 


n: Numbers of sequences of each group. 
*: Excluding Thermophis. 


Dipsadinae, and 7. baileyi vs. T. zhaoermii ranged 
from 0.229, 0.192 to 0.078 (Table 3). The number of T. 
zhaoermii vs. T. baileyi (0.078) was slightly more than the 
minimum number of Caenophidia (0.072). This lends to 
support the validity of recognizing T. zhaoermii as a new 
taxon (Huang et al., 2009). 


3.2 Topological results The topologies observed between 
the BI and ML trees were similar. Based on the ML 
and BI trees (Figure 1), Acrochordus and Achalinus are 
always subtended by the root node to the monophyletic 
Colubroidea. Viperidae formed a monophyletic group 
which was sister to a clade of Homalopsinae and Elapidae, 
as well as Natricinae, Dipsadinae and Colubrinae. The 
basic relationships of the commonly recognized families 
and/or subfamilies are similar to those of many earlier 
studies (e. g., Kelly et al., 2003; Burbrink and Lawson, 
2007; Wiens et al., 2008; Huang et al., 2009). 

Our data supported the monophylely of the dipsadines 
(99% in ML and 100% in BI) consisting of three 
subgroups: 1) “Clade Relicts”: American relict snakes and 
Thermophis (weakly supported here; Pinou et al., 2004); 
2) a strongly supported Central American group (92% in 
ML and 100% in BI); and 3) a strongly supported South 
American group (92% in ML and 100% in BI). South 
American and Central American clades huddled strongly 
and were supported by ML (74%) and BI (99%), and were 
sister to the root of “Clade Relicts”. 


4. Discussion 


One of the major challenges to inferring large-scale 
phylogenies is the expense and time necessary to obtain 
tissue samples and comparable character sampling for 
hundreds of species and many genes (Pyron et al., 2011). 
Under limited conditions, we only sequenced the ND2 
gene from two tissues (two species of Thermophis), 
retrieved 86 ND2 sequences from GenBank (74 
Caenophidia and 12 Henophidia), and combined them to 


G V (V/aligned sites) PI (PI/aligned sites) 
8.8 754 (73.1) 663 (64.2) 

8.9 525 (50.1) 319 (30.9) 

9.2 81 (7.9) 


Table 3 The uncorrected P-distances of ND2 gene, based on the 
data used in the present study. 


Taxa Uncorrected P-distance 
Maximum Minimum Mean 
Caenophidia™ n=73 0.356 0.072 0.229 
Dipsadinae™ n=15 0.255 0.123 0.192 
T. zhaoermii vs. T. baileyi n=3 0.078 


n: Numbers of sequences of each group. 

*, Excluding Thermophis. 

*; Including the sequences of Dipsadinae, Xenodontinae and North 
American relicts. 


infer a large phylogenetic tree. The relatively high support 
in main nodes (Figure 1) demonstrates that ND2 is a good 
genetic marker for inferring large-scale trees of snakes. 

Despite of the difference in nomenclature between the 
North, Middle, South American xenodontines and their 
subgroups by several researches (e. g., Pinou et al., 2004; 
Lawson et al., 2005; Hedges et al., 2009; Huang et al., 
2009; Zaher et al., 2009; Vidal et al., 2010; Pyron et al., 
2011), the relationships always cluster as: 1) xenodontines 
always consist of three parts, that is, Xenodontinae, 
Dipsadinae, and North American Relicts; 2) North 
American relicts served as outgroup of middle and south 
xenodontines; and 3) Thermophis was always placed 
together with (this study) or near (Pyron et al., 2011) the 
North American relicts. Huang et al. (2009) indicated 
that Thermophis is part of the subfamily Dipsadinae, one 
of the largest subfamilies of snakes with approximately 
90 genera and more than 700 species, with all current 
members of the subfamily are restricted to the New World. 
In the present study, the closest relatives of Thermophis 
appear to be the North American Relicts (Heterodon 
platirhinos, Farancia abacura). 

Some phylogenetic analyses discussed the origin, 
dispersal and faunal exchange of snake groups in the 
Holarctic Region (e. g., Vidal et al., 2000; Pinou ef al., 
2004; Burbrink and Lawson, 2007). Vidal et al. (2000) 
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Figure 1 The ML tree based on ND2 sequences of 87 species of Alethinophidia. The values above the line are bootstrap supports of ML, 
and the values below are posterior probabilities (PP) of the 50% majority-rule consensus tree from BI (%). Likelihood settings from best- 


fit model (GTR + I + G) selected by Akaike and Bayeslan information criterion in Modeltest 3.7 (-InL = 39597.2). Base = (0.4550, 0.3384, 
0.0405), Nst = 6, Rmat = (0.1017, 2.9511, 0.1697, 0.1903, 1.7759), Rates = gamma, Shape = 0.5219, Pinvar = 0.1577. 


No. 3 Xiaoyu SUN etal. 
inferred the phylogenetic relationships of the North, 
Middle and South American xenodontines from the 
sequence of portions of two mitochondrial genes (12S 
and 16S ribosomal RNA) in 85 species. They presented 
a hypothesis of an Asian—North American origin of 
dipsadines based on the monophyly of this group and 
their basal position with respect to North American taxa. 
In the present study, based on the analyses of the gene 
ND2, within the monophyletic clade of dipsadines (with 
supports of 99% in ML and 100% in BI), Thermophis 
clustered with the North American relicts to form a 
monophyletic subclade (despite with weak support), and 
rooted to Middle and South American xenodontines. 
This further suggested an Asian-North American 
origin of xenodontines. Similar topology (including 
Asian species of Plagiopholis, Pseudoxenodon, and/or 
Thermophis) was also revealed in the reconstructions of 
Lawson et al. (2005), Zaher et al. (2009), and Pyron et 
al. (2011). Extensive sampling (including Plagiopholis, 
Pseudoxenodon, Thermophis, and north, middle, south 
xenodontines) and same length genetic markers (including 
ND2) are needed in the future to obtain robust topologies 
with high support to further understand the history of this 
intercontinental group. 


Acknowledgements This research was funded by the 
National Natural Science Foundation of China (30870290, 
31071891) and the Students Science Research Program 
of Huangshan University (2010xdkj012). We thank Hao 
WANG and Fang WANG from Beijing University for 
helping in collecting samples. We also thank Baolin 
ZHANG, Chen LING and Jinmin CHEN from Huangshan 
University, and Tao WAN from the State Key Laboratory 
of Genetic Resources and Evolution, Graduate University 
of Chinese Academy of Sciences for their help in data 
processing. Likewise, we extend our thanks to Frank 
T. BURBRINK from College of Staten Island, City 
University of New York for his valuable comments and 
extensive revisions on this manuscript. 


References 


Alfaro M. E., Arnold S. J. 2001. Molecular systematics and 
evolution of Regina and the thamnophiine snakes. Mol 
Phylogenet Evol, 21: 408—423 

Burbrink F. T., Lawson R. 2007. How and when did Old World 
ratsnakes disperse into the New World? Mol Phylogenet Evol, 
43: 173-189 

Busack S. D., Lawson, R. 2008. Morphological, mitochondrial 
DNA and allozyme evolution in representative amphibians and 


Tibetan Plateau Relict Snakes of the Genus Thermophis... 167 


reptiles inhabiting each side of the Strait of Gibraltar. Biol J Linn 
Soc Lond, 94: 445-461 

Castoe T. A., Jiang Z. J., Gu W., Wang Z. O., Pollock D. D. 2008. 
Adaptive evolution and functional redesign of core metabolic 
proteins in snakes. PLoS ONE, Vol. 3 (5): e2201 

Castoe T. A., Gu W., de Koning A. P., Daza J. M., Jiang Z. J., 
Parkinson C. L., Pollock D. D. 2009a. Dynamic nucleotide 
mutation gradients and control region usage in squamate reptile 
mitochondrial genomes. Cytogenet Genome Res, 127(2-4): 112- 
127 

Castoe T. A., de Koning A. P., Kim H. M., Gu W., Noonan B. 
P., Naylor G., Jiang Z. J., Parkinson C. L., Pollock D. D. 
2009b. Evidence for an ancient adaptive episode of convergent 
molecular evolution. Proc Natl Acad Sci USA, 106(22): 8986- 
8991 

Dong S., Kumazawa Y. 2005. Complete mitochondrial DNA 
sequences of six snakes: Phylogenetic relationships and 
molecular evolution of genomic features. J Mol Evol, 61: 12-22 

Dorge T., Hofmann S., Wangdwei M., Duoje L., Solhoy T., Miehe 
G. 2007. The ecological specialist, Thermophis baileyi (Wall, 
1907) — new records, distribution and biogeographic conclusions. 
Br Herpetol Soc, 101: 8-12 

Douglas D., Janke A., Arnason U. 2006. A mitogenomic study on 
the phylogenetic position of snakes. Zool Sci, 35: 545-558 

Douglas D. A., Gower D. J. 2010. Snake mitochondrial genomes: 
Phylogenetic relationships andimplications of extended taxon 
sampling for interpretations of mitogenomic evolution. BMC 
Genom, 11: 14 

Guicking D., Lawson R., Joger U., Wink M. 2006. Evolution and 
phylogeny of the genus Natrix (Serpentes: Colubridae). Biol J 
Linn Soc Lond, 87: 127-143. 

Guo P., Liu S., Feng J., He M. 2008. The description of a new 
species of Thermophis (Serpentes: Colubridae). Sichuan J Zool, 
27: 321 

He M., Feng J., Liu S., Guo P., Zhao E. 2009. The phylogenetic 
position of Thermophis (Serpentes: Colubridae), an endemic 
snake from the Qinghai-Xizang Plateau, China. J Nat Hist, 43: 
479—488 

Hedges S. B., Couloux A., Vidal N. 2009. Molecular phylogeny, 
classification, and biogeography of West Indian racer snakes of 
the Tribe Alsophiini (Squamata, Dipsadidae, Xenodontinae). 
Zootaxa, 2067: 1-28 

Huang S., He S., Peng Z., Zhao K., Zhao E. 2007. Molecular 
phylogeography of endangered sharp-snouted pitviper 
(Deinagkistrodon acutus; Reptilia, Viperidae) in Mainland China. 
Mol Phylogenet Evol, 44: 942-952 

Huang S., Liu S., Guo P., Zhang Y., Zhao E. 2009. What are 
the closest relatives of the hot-spring snakes (Colubridae, 
Thermophis) the relict species endemic to the Tibetan Plateau? 
Mol Phylogenet Evol, 51: 438-446 

Huelsenbeck J. P., Ronquist F. 2001. MRBAYES: Bayesian 
inference of phylogenetic trees. Bioinformatics, 17: 754-755 

Jiang Z. J., Castoe T. A., Austin C. C. Burbrink F. T., Herron 
M. D., McGuire J. A., Parkinson C. L., Pollock D. D. 2007. 
Comparative mitochondrial genomics of snakes: Extraordinary 
substitution rate dynamics and functionality of the duplicate 
control region. BMC Evol Biol, 7: 123 

Kumazawa Y., Nishida M. 1995. Variations in mitochondrial tRNA 


168 Asian Herpetological Research 


Vol. 2 


gene organization of reptiles as phylogenetic markers. Mol Biol 
Evol, 12: 759-772 

Kumazawa Y., Ota H., Nishida M., Ozawa T. 1998. The complete 
nucleotide sequence of a snake (Dinodon semicarinatus) 
mitochondrial genome with two identical control regions. 
Genetics, 150: 313-329 

Kelly C. M. R., Barker N. P., Villet M. H. 2003. Phylogenetics 
of advanced snakes (Caenophidia) based on four mitochondrial 
genes. System Biol, 52: 439-459 

Ling C., Liu S. Y., Huang S., Burbrink F. T., Guo P., Sun Z. Y., 
Zhao J. 2010. Phylogenetic analyses reveal a unique species of 
Elaphe (Serpentes, Colubridae) new to science. Asian Herpetol 
Res, 1(2): 90-96 

Liu S., Zhao E. 2004. Discovery of Thermophis baileyi (Wall, 
1907), a snake endemic to Xizang AR, from Litang County, 
Sichuan, China. Sichuan J Zool, 23: 234-235 

Lawson R., Slowinski J. B., Crother B. I., Burbrink F. T. 2005. 
Phylogeny of the Colubroidea (Serpentes): New evidence from 
mitochondrial and nuclear genes. Mol Phy Evol, 37: 581-601 

Malnate E. V. 1953. The taxonomic status of the Tibetan colubrid 
snake Natrix baileyi. Copeia, 1953: 92-96 

Mulcahy D. G., Macey J. R. 2009. Vicariance and dispersal form 
a ring distribution in night snakes around the Gulf of California. 
Mol Phylogenet Evol, 53: 537-546 

Nagy Z. T., Lawson R., Joger U., Wink M. 2004. Molecular 
systematics of racers, whipsnakes and relatives (Reptilia: 
Colubridae) using mitochondrial and nuclear markers. J Zool 
Syst Evol Res, 42: 223-233 

Posada D., Crandall K. A. 1998. Modeltest: Testing the model of 
DNA substitution. Bioinformatics, 14: 817-818 

Pinou T., Vicarioa S., Marschner M., Caccone A. 2004. Relict 
snakes of North America and their relationships within 
Caenophidia, using likelihood-based Bayesian methods on 
mitochondrial sequences. Mol Phylogenet Evol, 32: 563-574 

Pyron R. A., Burbrink F. T. 2009. Neogene diversification and 
taxonomic stability in the snake tribe Lampropeltini (Serpentes: 
Colubridae). Mol Phylogenet Evol, 52: 524-529 

Pyron R. A., Burbrink F. T., Colli G. R., Montes A. N., Vitt L. J., 
Kuczynski C. A., Wiens J. J. 2011. The phylogeny of advanced 
snakes (Colubroidea), with discovery of a new subfamily 
and comparison of support methods for likelihood trees. Mol 
Phylogenet Evol, 58(2): 329-342 

Rao D. Q. 2000. Complimentary survey of the herpetofauna of 
Xizang Autonomous Region (Tibet) with discussion of their 
distribution and current status. Sichuan J. Zool. 19: 107-112 (In 
Chinese) 


Sambrook J., Fritsch E. F., Maniatis T. 1989. Molecular Cloning: 
A Laboratory Manual, Seconded. New York, America: Cold 
Spring Harbor Lab Press 

Thompson J. D., Gibson T. J., Plewniak F., Jeanmougin F., 
Higgins D. G. 1997. The CLUSTAL_X windows interface: 
Fexible strategies for multiple sequence alignment aided by 
quality analysis tools. Nucl Acids Res, 25: 4876-4882 

Townsend T., Larson A., Louis E., Macey J. R. 2004. 
Molecular phylogenetics of squamata: The position of snakes, 
amphisbaenians, and dibamids, and the root of the squamate tree. 
Syst Biol, 53: 735-757 

Tamura K., Dudley J., Nei M., Kumar S. 2007. MEGA4: 
Molecular Evolutionary Genetics Analysis (MEGA) software 
version 4.0, Mol Biol Evol, 24: 1596-1599 

Vidal N., Shannon G. K., Wong A., Hedges S. B. 2000. 
Phylogenetic relationships of xenodontine snakes inferred from 
12S and 16S ribosomal RNA sequences. Mol Phylogenet Evol, 
14(3): 389-402 

Vidal N., Maël D., David J. G. 2010. Dissecting the major 
American snake radiation: A molecular phylogeny of the 
Dipsadidae Bonaparte (Serpentes, Caenophidia). C R Biol, 333: 
48-55 

Wall F. 1907. Some new Asian snakes. J Bombay Nat Hist Soc, 17: 
612-618 

Wiens J. J., Kuczynski C. A., Smith S. A., Mulcahy D. G., Sites 
J. W. Jr., Townsend T. M., Reeder T. W. 2008. Branch lengths, 
support, and congruence: Testing the phylogenomic approach 
with 20 nuclear loci in snakes. Syst Biol, 57(3): 420—431 

Wüster W., Ferguson J. E., Quijada-Mascareñas J. A., Pook C. 
E., Salomão Mda G., Thorpe R. S. 2005. Tracing an invasion: 
Landbridges, refugia, and the phylogeography of the Neotropical 
rattlesnake (Serpentes: Viperidae: Crotalus durissus). Mol Ecol. 
14 (4): 1095-1108 

Yan J., Li H., Zhou K. 2008. Evolution of the mitochondrial 
genome in snakes: Gene rearrangements and phylogenetic 
relationships. BMC Genom, 9: 569 

Zaher H. 1999. Hemipenial morphology of the South American 
xenodontine snakes, with a proposal for a monophyletic 
Xenodontinae and a reappraisal of colubroid hemipenes. Bull 
Am Mus Nat Hist, 240: 1-168 

Zaher H., Grazziotin F., Cadle J., Murphy R., Moura-Leite 
J., Bonatto S. 2009. Molecular phylogeny of advanced snakes 
(Serpentes, Caenophidia) with an emphasis on South American 
Xenodontines: A revised classification and descriptions of new 
taxa. Pap Avulsos Zool (Sao Paulo), 49(11): 115-153 


